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ABSTRACT 


An algorithm for the application of a special purpose 
minicomputer to the automatic analysis and initial classi- 
fication of EW signals is developed. The minicomputer as- 
sembly language version of this algorithm is presented. A 
FORTRAN program which generates realistic data for signal 
analysis is described. Using the simulated data and an 
instruction level simulation program, the algorithm is 
tested using an IBM 360 computer. Microprogramming tech- 
niques are presented which optimize the minicomputer 


application. 
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I. STATEMENT OF THE PROBLEM 


Electronic Support Measures (ESM) are designed to 
extract information about hostile forces from the electro- 
magnetic environment they create. Since the earliest days 
of radar, there has been a parallel development of equipment, 
techniques and tactics designed to counter the use of radar 
by hostile forces by using the transmitted energy to detect, 
locate and track the hostile force, often without your own 
unit being detected by the searching radar system. 

The current proliferation of radar systems, telemetry 
systems, communications systems and navigation systems which 
all use electromagnetic radiation for operation has created 
a condition of overcrowding of the electromagnetic spectrum. 
Verceimportant from the ESM standpoint is that the détec- 
tion of a signal at almost any reasonable frequency is 
assured. What is not assured is being able to detect sig- 
nals which are determined to be of interest in a sufficiently 
timely manner to be of use to the Operational Commander. 

The dense signal environment which exists in world areas 
of interest, coupled with the rapidly decreasing allowable 
reaction times imposed by new and developing weapons 
systems, and the continuous requirement to maintain viable 
capabilities with a shrinking manpower force has imposed 
a requirement for faster, more accurate, automatic classi- 


fication of ESM signals. The development of techniques for 








their employment in the task of signal analysis in the ESM 
environment. 

Poudate.. the major effort has been to apply digital 
computers to the task of off-line signal analysis and 
specific threat identification problems. What is now 
required is the expansion of these or similar techniques to 
Momarestiye task Of routine surveillance over a wide fre- 
quency range, filtering received information for interest, 
determination of the parameters of received signals, 
initial classification and threat evaluation, and display 
Boman Operator Or Output to a remote location. All of 
this has to occur in real time with a high probability of 
eve lOpiie correct information. It is felt that such a 
system may be feasibie if it is designed taking advantage 
of the rapid data handling capability of digital computer 


systems specifically designed for this application. 





IT. CONCEPTUAL SYSTEM DESIGN 


A. SYSTEM OPERATIONAL REQUIREMENTS 

A system designed to fulfill the needs of automatic ESM 
surveillance classification and threat analysis has the 
following characteristics for effective performance. 

a. Broad frequency range. 

b. Simultaneous reception of all frequencies within 
the frequency range. 

eee dns Ot determining —rece1ved Signal characteris- 
tics within a reasonable accuracy. 

d. Ability to operate successfully in dense environ- 
ments with tens or hundreds of thousands of pulses per 
second being received. 

e. Ability to associate received pulses and other 
forms of radiation into pulse trains or other forms for 
identification and analysis. 

£. Means of determining the bearing of the received 
energy to within acceptable accuracy. 

Geeeerrecurave = meanseot Wdentifying =srenals as to Class 
and probable source. 

Mo oewtmes DrObaptt ity OL intercept. 


MmLGwemealsSe alarm race. 


Bs ole TEMSCOMPONENTS 
Conceptually such a system could be assembled using a 


series of Instantaneous Frequency Measurement (IFM) receivers 





that would cover the frequency bands of interest. The out- 
puts of these receivers would be applied to small special 
purpose digital computers which would perform the initial 
information filtering task of creating a data set of the 
parameters of all of the received pulse trains and exclude 
from further processing those signals which are obviously 
of no interest to the tactical environment. Information 
thus derived is then to be passed to a second computer 
System at considerably lower data rates which will allow 
detailed analysis of signals, final classification tasks 
and automatic threat evaluation. | 
Detailed design of an entire system such as that des- 


Splbea is beyond the scope of €his thesis. The problem 


the IFM and the associated digital computer and the infor- 


Haeronstiicerime task 1t 1s required to do. 


C. SPECIFIC SUBSYSTEM DESIGN ANALYSIS 

The specific design problem analyzed herein is stated 
as follows. 

Given an IFM recciver with an omnidirectional antenna 
and reasonable sensitivity and probability of detection at 
ailerrequencies within its band width, analyze the feasi- 
Dele sot PrTOvidine a special purpose digital computer that 
will accept all information developed by the IFM and provide 
as real-time outputs the frequency, pulse width, pulse 
repetition interval, scan type ard-the time of last 


Imcercept. 
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II1. COMPUTER SELECTION FOR PURPOSES 
OR Pens retin PY ANALYSIS 
To complete the analysis of the feasibility of the auto- 
matic IFM system, it was desirable to actually encode and 
Simulate the operation of the system using an existing com- 
puter instruction set. It was also of interest to determine 
the required execution times, pinpoint improvements to the 
basic algorithm, establish required additional ‘special 
case’ instructions which could improve the efficiency of 
execution, and establish desirable hardware arrangements 
Which would enhance the overall operation of the computer. 
Through the experience of operating a simulation program 
and the test of various software and hardware techniques, 
the detailed specifications for a special purpose computer 
can be determined. This computer would be specifically 
designed and optimized for the demanding high data rates 
anticipated. Investigation disclosed the existence of a 
small high speed special purpose computer produced by 
Applied Technology Corporation (ATI) in Sunnyvale, Califor- 
Mia. This computer is designated as the Applied Technology 
Airborne Computer (ATAC) and is a 16 bit fixed point 2's 
complement machine. A number of unique features of this 
computer made it an obvious choice for this study. The 
general computer characteristics and these special features 
are summarized below. Further detailed information is 
available in ATAC, Applied Technology Airborne Computer, 


eet. Principles of Operation. 


itt 





1. 16-bit instruction word length with a 437 nanosecond 
mejor cycle time. 

Z. Sixteen general registers. 

3. Four basic instruction addressing modes with three 
supplemental modes allowimne a wide range of direct; indexed 
and immediate execution instructions. 

4. Asynchronous memory with provisions of up to 4 
Direct Memory Access channels. 

5. Memory types can be mixed in any combination of Read 
Only Memories (ROM), Bipolar Ramdom Access Memories (RAM), 
Read Mostly Memories (RMM) and slower core memory devices. 

Sa lwelVe Maan prior cy) sneer iupts. 

7. “lO>-bit = paraliel Input cutput operations. 

Sa EXtremely isitalh sa gzemonGge pOWemmuculi hememecian© 
pounds, 115 cubic inches, and a power consumption of 
approximately 70 watts for a nominal configuration. 

9. An existing software support package designated as 
the AITAC Programming Support System (APSS). 

10. Microprogrammed instruction set allowing ease of 
modification of the instruction set within the hardware 
Minitacrons of the machine. 

11. Computer assisted hardware design and fabrication 
allowing relative ease and speed of modification of the 
Hoteware Contiguration (within Jamits) to permit optimiza- 
iGmoneor the Computer for specifire applications. 

Mines features listed und@y 9, 210) and Tl are the essential 


features required for this feasibility study. As the basic 


I 





speed, memory configuration and programming flexibility 


required are provided by the existing ATAC, arrangements 


were made with ATI to allow use of portions of the APSS for 


this effort. 


A. ATAC PROGRAMMING SUPPORT PACKAGE (APSS). 


APSS is a computer simulation program designed and pro- 


grammed for use on IBM 360/370 series Operating Systems. 


The complete package consists of the following items. 


Volume 1 - 


Volume 2 - 


ATAC BASIC PRINCIPLES OF OPERATION 

This document describes the basic ATAC machine 
CharacteristicSsana! ehe Mrstruetion set. 

ATAC PROGRAMMING SUPPORT SYSTEM (APSS) 

APSS is the support software system for ATAC 
which allows preparation and debugging of 
programs via an IBM 360/370 computer system. 
The dOCUMeh ES MI Chis eveliMesg@esen Tic nO LO 
set up and use the APSS system and includes 
publications describing the PL/ATAC Compiler, 


Assembler, Loader Simulator and microprogramming 


faecal t Lece 

SeGt Ont - AFPSS Job Meanacement 

Section II - PL/ATAC Compiler 

Section 111 - APSS Assembler/Loader 

Sect zoned Vv = HPoom Stile rcion evel Simulator 
Secrlonn - APSS Microprogramming Software 


ILS 





Volume 3 - ATAC SYSTEM DIAGNOSTIC PROGRAMS 
This document summarizes the use of ATAC's 
Diagnostic programs for instruction test, I/0 
test, Imterrupt systenmeheckout, and condition 
code. test. 

Volume 4 - ATAC CIRCUIT BOARD TEST PROGRAMS 
This document summarizes the ATAC circuit board 
Pec proceaunese 

Vemume S = ATAC OPERATOR'S GUIDE 
Procedures “orm ence Opera eronwot an ATAC 
computer are described in this document. 

Volume 6 - ATAC 16A MICROPROGRAMMING AND SIMULATION 
PROoGeduGes tOmedcey> lopiicmleroproct alc mse, 


the ATAC 116A Central processing unit. 


PicmoOorulons Of APSSsUtilizedqmropellic study are dos ted 
betow. Other volumes were not required but their existence 
as noted. Future efforts in the development of specifica- 
tions for a special purpose computer should be made with 
the knowledge that such systems exist. The advantages 
accrued by having a computing system specifically tailored 
to the task at hand are self evident. Additional flexi- 
bility 1s provided through the use of microprogramming in 
that a completely different instruction set can be defined 
by replacing the microprogram Read Only Memory allowing 
the same basic machine to be optimized for several different 
tasks and reconfigured by changing a small inexpensive 


module. 
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Volume 1 - ATAC BASIC PRINCIPLES OF OPERATION 
Volume 2 - ATAC PROGRAMMING SUPPORT SYSTEM (APSS) 
Secrr1om | APSS Job Management 
Section III APSS Assembler/Loader 
section IV APSS Instruction Level Simulator 
The APSS Instruction Level Simulation Program is 
currently a resident program on the IBM 360/67 at the Naval 
Postgraduate School. Details concerning the implementation 
of the program are contained in Appendix A. Details con- 


cCérning its use are contained in Appendix B. 


ae: 





IV. THE ITFM DIGITAL PROCESSOR 


The initial design concept for the interface between the 
IFM and the special purpose digital computer is a common 
memory storage area to which both units have Direct Memory 
Access on a non-interrupt basis. The data to be passed 
between the IFM and the computer consists of the frequency, 
pulse width and time of arrival of all received signals. 
Using this information, the computer can classify signals 
as pulsed type radar signals based on the indicated pulse 
width. Further analysis is conducted by establishing 
initial pulse files segregated by frequency and ordered by 
the sequence in which they were received. An estimate of 
the PRI is made by averaging the time between the arrival 
of adjacent pulses in the file. Small variations in the 
PRI estimate are treated as noise. Large variations in the 
time between pulses are considered caused by antenna scan 
@iaraecteristics. Variations between these extremes are 
fomsotdered due to the superposition of multiple pulse 
Geis Of due to a complex PRI “yatter' imposed on the 
emitter. Variations due to noise and antenna scan are 
Pwateamrncernallys Variations due to interleaved pulse 
trains and complex jitter patterns are signalled and further 
analysis of these cases is handled by the main classifica- 
tion computer. The design calls for the special purpose 
IFM computer to output the following information for each 
Meese train received. 


16 





a. Frequency 

b. Pulse width within limits, the value 

c. Pulse width out of limits, a designation 

d. Stable PRI, the value 

e. Unstable PRI, a designation 

f. Antenna scan time 
Time of last update 

hy. “Sinrenmordered »wampiiise train®of speeiwiied enprh to 
be used for analysis of complex signals by the main 


elassitication computer. 


A. TIME ENCODING 

Detcrmination=oe: thee pullse wadth and the tine of arrival 
Pedqmires 4a digital clocking system whose output is directly 
translatable into real time. For purposes of analysis it 1s 
assumed that when the IFM threshhold is exceeded in a 
specific frequency band, the current value of the clock is 
loaded into the time of arrival words. Simultaneously the 
pulse width word is cleared and then incremented at the 
mockerate. Whem the pulse Energy ceases, the pulse width 
Vetmmeomltains =the pulse width scaled by the clock rate. 

The clock consists of a 40 bit counter. The clock rate 
fi UsMerahertiz. This corresponds to a least significant 
' bit value for pulse width of 0.1 microseconds. The pulse 
Micemevordeconsists of bits 0-15 of “the counter which 
allows representation of pulse widths of up to 6.5 milli- 


seconds tre@ting the most significant bit as a magnitude 


ee 





instead of a sign bit. The time of arrival double precision 
word consists of bits 8-39. This corresponds to a least 
Significant bit value of 25.6 microseconds. The maximum 
cock time (time to overflow) is 15.3 hours. The most 
Significant bit of the clock is treated as a sign bit to 
allow detection of the data out of sequence condition. 
Accuracy is determined by the stability of the basic 
Clock oscillator. Measurement is made to within 0.05 micro- 
seconds for pulse width and to within 12.8 microseconds 


for PRI. 


B. THE COMMON MEMORY BLOCK 

The common memory block is used as a large buffer to 
Matenm the differing data rate handling capability of the 
Bimand the computer. For purposes cf analysis, it ls 
assumed that the IFM can respond and prepare a series of 
four data words in the times achievable using TTL logic 
G@mepeuatry. this allows the IFM to load the double pre- 
cision words for time of arrival and the word for frequency 
first, the word for pulse width as soon as the pulse has 
completely been received and accomplish the entire transfer 
of the four words in approximately 3 microseconds or the 
pulse width of the received signal + 0.5 microseconds which 
en molomeer. For signals with pulsSe widths of less than 
3 microseconds, the IFM data will have a maximum burst rate 
Sieapproximately 300 KHZ. 

Simulations conducted using the algorithm described in 


Peeoer IY indicate the averagé time required for the computer 
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to process each received pulse including read and write 
times 1S approximately 45 microseconds with improvement to 
35-40 microseconds considered readily possible. This 
corresponds to a maximum average data rate of approximately 
22,000 pulses per second. 

The IFM loads data words into the common memory block 
sequentially starting from the '0' location and wraps around 
Vieimene block 1s full. Also starting from the '0' location, 
the computer accesses the data words and associates them with 
Peep UGO pilates pulse traineand then accesses the mext aord. 
The size of the buffer establishes the length of the time 
that transitory high data rates can be sustained by the 
‘rem before the conputer 1S ‘lapped' by the IFM input. 

In this case, the computer is always accessing the most 
Becent data set and the information stored by the IFM pre- 
viously is lost. When the data rate decreases, it may 
Peeores possible tor the computer to catch up to the I&M to 
Pemove the lap condition. In this case, the data read by 
the computer will be out of sequence and will be discarded. 

Pulse information is stored in the common memory block 
iMmmamscwiles Of Lour words per pulse. These words are 
defined as follows. 

iarcaleney. Frequency 15 represented by an integer 
fOrresmoemamme to the frequency bin in which the detection 
OCcurs 

beer nlse width. The pulse width is represented by an 


Meeaegrewhiich is a scaled time value. 
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c. Time of arrival. Time of arrival is expressed as a 
fomere —preetston integer, the lower half of which is the 
third word, the upper half the fourth word. | 

The computer can keep track of the existence of new 
data in the word block associated with an individual pulse 
by resetting the frequency word to 0 after it is read. By 
testing the frequency word, the computer establishes whether 
the IFM has loaded new data since its last look at this 
location. 

By making the IFM output and the computer input com- 
pletely independent and asynchronous in operation, timing 
and interrupt communications, with their high time overhead 


Posts are avoided. 


C. THE IFM COMPUTER ALGORITHM 

An algorithm has been developed which defines the oper- 
eeromoerequited to develop pulse trains and extract PRI 
mieormatzon trom the [FM output data. This algorithm 
causes the computer to read a memory location and determine 
leader 1S available. If the data is available, the fre- 
quency word is used as a relative address to enter a file 
which provides linking information for the pulse train 
mibeomacvetopea Lor that frequency. This file data word is 
utilized as two half words. The lower half is used as a 
Dulse train file pointer. the upper half is used as a 
Polilecre1Or an Overtlow storage area in the event there are 


two Or more pulse trains associated with a given frequency. 
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The primary means of pulse train association is by frequency 
discrimination. 

The pulse train file pointer is used as a relative 
address for entry ito amcontactmiilesentry tablegws. The 
Giplewentry Corresponding to the pulse train file pointer 
contains two words. The first word contains the frequency 
Siecie pulse train for cross reference purposes. The second 
word contains the absolute entry address of the actual file. 

The pulse train file is used both to accumulate data for 
calculation of parameters and to provide output information 
tOoetne main classification computer by using Direct Memory 
Access techniques. The file is a block of 64 memory loca- 
tions which are treated as 16 four word groups. This allows 
Sroneme 2ntcrmation on 16 consecutive pulses which are used 
to develop an estimate of the PRI. Before 16 pulses have 
been received, the file is said to be building. After 16 
pulses have been received and the PRI estimated, the file 
memcaed teoebe closed. If no current pulse information is 
eomieeined in the file, 1f 18 said to be empty. The distinc- 
Prone ts requiread because the contents of the words of the 
file have different meanings for the different cases. The 
Cemmemes OL the fille words are described below. The choice 
of 16 pulses used for PRI determination is arbitrary and 
additional testing and evaluation is required to determine 
Promo pulmunm number of pulses considering detection proba- 
bilities, normal antenna scan rates, false alarm rates, and 


Sii@eatron and operational experiences. 
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IT rile Deseidang . 
a, word 1. SPulse width of the initial pulse 


Wewectca When the pulse train file was established. 

be Word=2. Colmeer imareatine the number o1 
pulses Storen anernis file. The counter counts from -16 
wor 0). 

CmENOGG. 5,64 es Deuble precision tame Gof arraval 
PEseiemrirst pulse. 

aememicrd 5. Pulse width of the second pulse. 

e. Word 6. Time difference between the arrival of 
the second pulse and the first. 

iW Org =), of Doublew precision Crime Ob arrival’ of 
emer second pulse. 

ge. oucceeding groups of. 4 words are arranged ln 2 
Similar fashion with the second word of the group containing 
tiem eame Interval between the pulse and the one preceding 
em 

oe enc LOSed, 

aumeniomd 1. Average pulse width of the file. 

beg Woma 2. “positive integer represents the 
average PRI: If the location contains a zero, a stable PRI 
Camot be Ggetermined. 

e- Nord 5, 4. Double precision time of ast imtror- 


Macion update. 


Ze 





d aeilior dus, 6. Double precaisaien time of armival of 
the first pulse of the latest consecutive group of pulses to 
Bbesreceived, 

e. Word 6, 7. Double precision antenna scan rate. 

Cre elemennty.: 

The file empty condition is indicated by the first 
(meoonrne COneact Lave entry file being zero. The con- 
Belles OL the pulse train file Fhoulemoe considered as being 
unspecified. 

The scheme of using the frequency as a relative 
Steuyeaaaress to the frequency fle, Wsitg the data comtammed 
in the frequency file as an entry into a contact file and 
Doneeehesdata contained in the contact file to prowide the 
Mela waddiress of the pulse train files is required tommein- 
tain pulse train files without the requirement to have a 
separate pulse train file for every possible combination of 
frequency and pulse width. A total of 64 pulse train files 
can be established and maintained simultaneously in the 
Current configuration. 

In the event that two pulse trains of differing pulse 
widths are received in the same frequency bin, provision 
has been made to discriminate these signals by pulse width 
and establish and maintain separate files on these signals. 
If interleaved pulse trains of nearly equal pulse widths 
are received on the same frequency, the algorithm will con- 
sider the pulse train as having a complex PRI and enter a 


Senonin tre PRI word: The main classification computer can 


Li 





themmaesignatemthe requimement to establish a #ile of 

Supe lecnt@ length eto colfecteiough pulses for®analysiS by 
inserting the required number of pulses into the PRI word. 
After filling 60 of the 64 assigned locations for the pulse 
train file, the remainder of the long file is inserted into 
an overflow section of memory. Words 61-64 are used to pro- 
Vide wimtiinrg = 1nsoriat lon . 

Calculation of the PRI after 16 pulses have been received 
1s accomplished by averaging the time interval between pul- 
ses. Sixteen pulses are contained in the pulse file and the 
seventeenth pulseto arrive causes the start of this calcu- 
lation to provide™¥6 time intervals for the averaging. Tie 
mieervyals are treated in Teverscrcrdeyeao. thelr apperyanec. 
The latest Interval is treated as a standard and the pre- 
vious interval is compared with it. If the two intervals 
dmemiot Medrly Cdualy= the conditrommen one being nearly 
equal to a multiple of the other is check for multiples up 
to 4. If this condition is met, both intervals are made 
equal to the smaller interval and their sum is accumulated. 
PhemceiecmecOnGltron 15 NOu met, bOeMwInmtervalis are still made 
equal to the smaller interval (the longer interval discarded) 
their sum accumulated and a discard counter is incremented. 
Ditemprocess 15 repeated with the other time of arrival 
increments in the file until all 16 intervals have been 
summed or until the discard counter reaches 4. If 4 pulses 
are discarded, the PRI is considered unstable and a zero is 


ansertcd into the PRI word. If ali 16 pulse intervals are 
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summed, the resulting sum is divided by 16 (shifted right 
4 bits) and the result inserted into the PRI word. Con- 
currently, the pulse widths are also summed to calculate 
the average pulse width of the file in a similar manner. 
After the file has been closed, succeeding pulses 
received are Uti lizceeromuUndarestnecedacasim: the file. “lta 
received pulse width is within limits of the file pulse 
width, it is compared with the file pulse width to determine 
its relative size. If the new one is smaller, a1 is sub- 
tracted fromthe file pulse rene li lagoceraeebeis 
added. lIlhis featurevensuresmtnat lone termedrrcts on tne 
parameter value can be tracked. If a PRI has been established 
for the file, and the time interval between succeeding pulses 
exGeeds Cscablished limits Min omncwepilse 15 Considered to 
have established simple antenna scan parameters and an 
anetina SCan rate Cote uomcalcmlatcd, 
The IFM computer algorithm is shown in flow chart form 


in Figure 1. 
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V.. DATA SIMULATION PROGRAM (DSP) 


The DSP was written by the author to provide a relatively 
realistic input data environment to the ATAC Simulater. The 
program was written well before the ATAC Simulater became 
available and includes many features which are not useful in 
in eet EReNntoannolication. 

(Hees spreovr~des Simulated signal levels arriving at an 
IFM receiver in a dynamic environment. Input data provides 
Mitaal set up positions for the IFM and up to 10 emitter 
Pl aeeorms. The IFM and all of the 10 platforms are given 
courses and speeds which are used to update their positions 
as the simulation proceeds. Each platform may have up to 5 
emitters on board of any specified frequency, PRI, PW, 
anetnna scan, antenna gain, and antenna beam width. The 
program establishes the initial antenna position by using a 
uniformly distributed random variable, updates the antenna 
position, calculates the Sipnal strength of each transmitted 
pulbceswien received by the [FM and outputs those pulses 
which are above the specified threshold sensitivity of the 
Lim lre OFUuEput Consists of the time sequence of pulses 
Feceived by the [FM including the frequency, pulse width, 
time of arrival, bearing and signal strength. 

The data used by the ATAC Simulater consists of four 
Mireeers £or Gach puls@. These words are: 

Im Frequency bin. The bins ar€ numbered in ascending 


Mecer sor frequency from 1-128 although the only restriction 
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is that there has to be less than a total of 255 bins for 
the current ATAC Simulater program andthe 0 bin cannot be 
used as it is used by ATAC to signal a no data state. 

2. @Bulse width. The pulse width word is an integer in 
the range 0-65,535. Any time scaling factor may be used. 
The current convention adopted is that the integer 1 cor- 
responds to a time interval of 0.1 microseconds. This 
alikows representation of any pulse wadth from 0-6.5 mzlii- 
seconds to within 0.05 microseconds. 

3. Time of arrival. Time of arrival is expressed as a 
iouple precistom (S2ebat) integers Worass) 1s. the Least 
Significant half of this enteger and holds time values up 
tom0:.84 seconds with the dsea@st significant bit comresponding 
ieme2 a. O MLCrOSeconds. 

4. Time of arrival. The upper half of the time word 
holds time values up to 15.5 hours. This clock update 
rate and time measurement accuracy is considered adequate 
for contact analysis and PRI measurements. 

The DSP generates data in the above format at a rate of 


approximately 200 pulses per second. 


Doe ene@mPLES OF OPERATION 

1. Position information. All position information is 
initially supplied in X-Y grid coordinates in nautical miles. 
Velocities are input in knots and direction of travel in 
degrees clockwise from North. The DSP converts this infor- 


Meron to X-Y coordinates in yards réllative to the IFM and 
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relative velocity components ay the X and Medirections in 
yards per second. The Relat ivceVelocityv,. Closest Pommteot 
Peproach (CPA) and distance to CPA is calculated for all 
platforms. 

CoeeEMittcumanpue data consists of the transmit frequency 
fd2. pulse width in microseconds, PRI in milliseconds, 
transmitted power in KW, antenna gain in db, antenna beam 
width in degrees and antenna scan rate in seconds per revo- 
iition. DSP converts ald input units to common units for 
computation. 

3. IFM input data consists of the maximum and minimum 
frequencies of its pass band in MHZ, the minimum measurement 
time of the IFM in microseconds and the threshhold sensiti- 
ioe db. 

tewetilazes the input data to calculate the received 
Signal strength for each emitter as a function of range and 
the angle between the emitters antenna and the IFM platform. 
The maximum range of detection assuming maximum antenna gain 
is calculated and compared with the initial range and the 
CPA to determine when the emitter will come within range of 
the IFM, if ever. The pulse repetition cycle and antenna 
initial angle are assigned by assuming a uniform distribu- 
Bieoeesbarting positions in these cy¢lles and accessing a 
Umitorm distribution random number generator. When the 
Sieeeredis within range, the portions of the antenna scan 
cycle that the antenna gain is high enough to allow detec- 


mmon are calculated. 
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On the basis of the above calculations, DSP stores the 
time of arrival at the IFM, counting transit times, for the 
first pulse of each emitter which will be received with a 
Signal strength equal to or greater than the receiver 
sensitivity. These times are stored in an array which is 
then searched to find the next occurring time of arrival. 
As each pulse time of arrival becomes the next in the 
sequence of the receiver pulse train, DSP calculates the 
time of arrival of the next pulse from that emitter and 
stores that time in the array. Each time update continues 
to consider the transit antenna scan and the relative 
motion between the IFM and the emitter platform. 

The DSP currently models the IFM with a 100% probability 
Bt aetection of all puises arriving with signal strengths 
equal to or greater than the threshhold sensitivity. 
Additional coding can readily be added to provide a less 
than perfect detection scheme, false alarms, and an IFM 
mirectional antenna. 

The DSP has been provided with selection indicators for 
selecting random variations to all emitter parameters. Cod- 
ing has not been included to accomplish any such variation 
But Can be readily added. 

Tie USPecurrently imciudes only two options for the 
emitter antenna scan type--steady and circular scans. Cod- 
ing for additional scans can also be readily added. 

The overall flow diagram for the DSP is shown in 


mugure 2. ¢ 
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Data Simulation Program Flow Chart 
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wel, )|©6 «COP TEM DOALTION OPSiHE PPM DPGPPAL PROCESSOR 


The ATAC Assembly Language program listed in Appendix C 
and described in detail in Appendix D is written using the 
Supplied instruction set. The major advantage of a micro- 
programmed computer is that the instruction set can be 
plored for weprimunme@ertie@iency inWexecution of a specific 
algorithm. Two primary results of this optimization are the 
saving of instruction fetch times when a series of instruc- 
tions are microcoded and the time savings realized by 
performing operations simultaneously when the data paths do 
not overlap. Modifications to the instruction set are 
constrained only by the hardware configuration of the 
computer and the programmer's imagination. 

In addition to the time savings associated with modify- 
ing the instruction set additional benefits can be accrued 
by modifying the hardware to allow simultaneous input and 
output operations, buffering the memory interfaces with 
heel speed registers and instruction intemleaving techniques. 

Modifications of both microcode and computer hardware 
can be tested and evaluated by using the APSS Microprogram- 
ming Software. By simulating the desired microcode sequences, 
debugging and optimization efforts can be expended before a 


commitment to actual design change is made. 
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A. MICROCODE SEQUENCES 

Serres Of ianstructions which are shown to use a large 
percentage of the execution time of the program should be 
examined for microcoding first. The time expended by each 
instruction is output by the APSS Instruction Level Simula- 
tor in the form of a time histogram. An example of this 
histogram is shown in Appendix D. This histogram shows that 
the two instructions which use the most time during the 
execution of the program are the store data and load new 
data instructions. The program sequence which ligalsbotrstexs 
these instructions is required for every input pulse and 


ms listed below. 


Prstruictaon Purpose ee Qinet Omer dic 
(microseconds) 

i DRS 4 Store data Des 

BbRGo-U, TOP Branch to beginning Mey ol) 

ee 15. 054 Increment address Bes 

Rel DR,O,4 Piet ee WeGia td Chas. 

EVR -RyS, 1 Move frequency word BT! BIS. 

BRGor cL LOr Feserromipteseiee OL «data 2457 5 

ILJDNe ss a0 Pipe nie ners cer, | ue es 


oR ekxXs 1.051 Score was) wnethescata DLOcIA T5025 
AOD 1S. 0.4 Reset address pointer BL igs 
Oke -DAcO CRESS Fetch frequency file entry 1.3125 
PES IZ. SUVARI Test tor aspulse train file .4375 
Total “tine sexoenaed eee a2 
The average time expended per input data set during 


Sammy iation runs was 45 microseconds. It is seen that 303 
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Paeenisetime 1s 8spent in the store/fetch phase. This 
sequence placed in microcode would save an initial 4.3 
microseconds in instruction fetch times. Additional 
Savings are possible by the simultaneous incrementing of 
Ene waddress registers agd simultaneous input and output 
operations as the input data and output data are using 
independent segments of memory. It is estimated that the 
SxCCURM@MEIMe Of Mihe Migmocoded sequemce woulla be appmex- 
mately 6G smicroseconds which reduces the program execution 
mIncaMen Pulse to appreximatedy 40 smnierosecondss = Addational 
microcode sequences can be developed, especially for the 
routine update data section of the program. A time saving 
Pieovlicn 1s probable. | 

Using the technique of microcoding instruction sequences 
and Optimizing the microprogram data flow, the average time 
PepeapulicemGcame pe reduced GO co1tO 50 Microscconds. Phas 
Step alone increases the average pulse density capability 


Hianomemaystenl tO o9h -4°0Kepulses per second, 


B. HARDWARE MODIFICATION 

tie time required to access ama “store the data can be 
EeiueciwupDy a tactoOr Of Lour by altering the means of storing 
tt.) lt the data is stored in four independent memory 
nMeGwbescnmparaitel transfer of the four words is possible. 
A hardware modification is required to allow the computer 
PeOmmomamueur resisters in parallel by supplying the addi- 


Priimeudata paths required, 


S5 





Supplying an input data buffer and an output data buffer 
each of four words would allow loading the input buffer 
during execution of instruction sequences in the middle of 
the loop and then loading the working registers in parallel 
im one major cycle. For storing data, the buffer can be 
loaded in parallel in one major cycle and the data stored 
in memory from the buffer during execution of the remainder 
of the program. If two buffers are provided both the fetch 
and store operations can be conducted simultancously. From 
the program's point of view, the execution time required for 
the 13.8 microsecond sequence of V.A above becomes .4375 


microseconds. 
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Ve aa ieee Ne oro Or bade I FM yiGrraL PROCESSOR 


Computer simulation of the IFM digital processor has 
Mecheconauerea.. | !ne=aleorithmn presented in Chapter III was 
coded in ATAC Assembly language. The resultant program is 
described in Appendices C and D. The Data Simulation Pro- 
gram (DSP) was used to provide a realistic environment of 
received pulses in order to test the concepts presented in 
Chapters II and III. Sufficient time was not available 
to conduct simulation runs to test all conditions and to 
develop more than an initial evaluation of the concept's 
Viability. A summary of results obtained, conclusions and 


recommendations for further study are included below. 


A. SUMMARY OF RESULTS 

The ATAC computer was modelled as provided for in the 
Aros Imstruction Level Samulator usingw@the program de- 
scribed in detail in Appendix C. Due to limitations imposed 
iy wne Simulator, data was inserted into the ATAC construc- 
MiVemnhelor wepr lor toluene DecInning Otmcamulataon. The data 
was prepared by the Data Simulation Program described in 
detail in Appendix D. A total of 1024 radar pulses were 
utilized in order to check on the accuracy of pulse train 
building, the development of PRI information and the average 
EMicetcduared ste nandle cach pulse. The results of the 
Simulation were the completely accurate discrimination of 


pulses by frequency and pulse width and the calculation of 
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mie correct PRI in all cases. The total time required for 
the ATAC to process 1024 pulses was 0.046 seconds or an 
average of 45 microseconds per pulse. The results are 


summarized in the table below. 


Simulation input data. 


Number of radar emitters 9 
io Gale iuinne? Omi LSeseitea liaed 1024 
Total real time simulated 0.1568 sec 


eamulation results. 


Accuracy of pulse association 100% 

Total ATAC simulated computer time O20 4 Gc ec 
Average Lime per pulse as M1 €rosec 
Maximum average data rate possible 22,000 pps 


fe =CONCLUSIONS 

The simulation was conducted with stable pulse widths 
and PRI. Even so, on the basis of the results of the simu- 
feati1on and on the information contained inChapter V, it is 
feasible to construct a system capable of significant amounts 
me alcomatic Signal classification in a high density environ- 
Ment With the average rate of intormation in the bandwidth 


of an IFM in the range of 25K to 50K pulses per second. 


~C. RECOMMENDATIONS 
Based on the results of this study, the following reconm- 


mendations for further study are submitted. 
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1. Continue simulating the high density environment with 
the programs listed in Appendices C and D adding to the 
Somplexity of the input data in order to determine the opti- 
mum algorithm for the solution of the problem. 

z. Use the results of the simulations to develop 
additional special purpose instructions which optimize exe- 
cution of the programmed algorithm. 

3. Use the results of the simulations to develop 
désirable hardware reconfigurations which will optimize the 
execution of the programmed algorithm. 

4. Use the features provided in APSS to test and evalu- 
ate the newly developed hardware configuration and instruc- 
mron set. 

b.  Expand@™the Wasic preecedures Gutlenred=nerein™@tro other 
melas that encounter Wigidata “rates and require automatic 
data processing. Examples are radar systems, communications 
systems and sonar systems. The desirability of specifically 
designed processing equipment has been demonstrated. 

Surrent technology is producing these processors which 
prilinate the processing time and price penalties paid when 


MeE1lizing general purpose digital computers. 
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APPENDIX A 


INSTALLATION OF THE ATAC PROGRAMMING SUPPORT SYSTEM (APSS) 
AT THE NAVAL POSTGRADUATE SCHOOL 


meine ATAC Programming Support System operates on a ecail 
or memote IBM 360/370... 
ATAC Principles of Operation (Volume 1) 
=lfeysomething can possibly go wrong, it will... 
Murphy's First Law of Dynamics 
A. Installation of the ATAC Programming Support System on 
the IBM 360/67 Operating System at the Naval Postgraduate 
menool by the author was Gece d Ilwerdear to allow pRo- 
gramming of algorithms in an actual representative machine 
assembly language and to determine execution times and over- 
load points for the algorithm/computer team in an EW 
mivironment. After negotiations with AT] Corporation, 
arrangements were made for a tape copy of the portions cf 
Poomuecdilyed ito, conducteanstruction levelesinulations. 
Smlicw cape was received on 109May 1974. The ensuing events 
Mecmrpmccented aS an example of interfacing problems that 
mie —CHCOUnLeErCa. sProblems encountered during the program 
mipkementation are listed in the ordéergef their occurrance 
Peon? With the correcti ves aqewuon Gakens 
ea piemdens ity . 
a. dhe original tape was prepared by ATI using a 


remote terminal to an IBM 370 installation. This tape was 
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taken to NPS by the author and an attempt was made to copy 
pronetice tape. ine tape Simply ran out indicating that the 
Operating System was unable to find any of the specified 
files On the tape -aeeutility progwameeAPEOUT’ waise used to 
determine what, if anything was on the tape. According to 
the utility program output there was nothing understandable 
necoOLded Olwtne tape at all Agbelated anvestigarrvon of the 
JEL -iInStruccioOns that prodauced™ tne tape revealed that it 

had been prepared on a dual density tape unit and recorded 
with a density of 1600 baits per anch. The installation at 
NPS does not support recording densities of greater than 800 
|pjokn ie 

b. Corrective action. A second tape was generated 
by ATI specifically recorded at 800 bpi. This tape was 
picked up by the author on 14 May 1974. 

2. Block length. 

a. Files 1 and 2 of the new tape were successfully 
read by the NPS computer. These files contained instruc- 
tions on how to access File 3 which contained the Instruction 
Peover oimulator mr object™module form. The required steps 
Wel cr 

1) Copy the object modules from the tape onto 
a direct access storage device such as a disk pack. 

2) Build an absolute load module by providing 
local linking data between the various portions of the pro- 
Maan COmtained in the Object modules. In object module 


One tre rogram COnSists OL a group of subroutines which 
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are stored with instruction and data locations specified 
as relative to the individual subroutines start point. 
Building a load module orders these individual sub programs 
and establishes a common relative base address so that when 
mars Cntry point is specified (when loading the program for 
Eee eCUtION) the instruction sequence proceeds as omiginally 
designed, including linking between these sub programs. 

3) Load and execute the APSS Simulator. 

Step 1 was accomplished. However, when attempting to 
meacecas the file created by copyam@ethe tape, the computer 
feeminanple to @taentify the contents of the file and returned 
fic error message ‘Permanent 1/0 Error on Attempting to read 
Prile Headings'.” Investigation indicated that this could be 
mrsca jy an incompatmple™ block™ieneth of the original 
object modules. Due to size limitations and buffering con- 
ractations, records (data) transferrec by computers are 
Oeeanized into blocks of data. These blocks are iden#@ified 
by length and where the next sequential block of data is to 
[eround. lo Make sense to a computer, when the data con- 
Sists of an executable program, it must be reassembled with 
mee same block structure in the proper alignment that it had 
when the program was assembled. Investigation revealed that 
the original assembly of the program was conducted using a 
Soo disk pack and was creatéd with a block length of 13k 
mes, ihe maximum block length supportable at NPS is 7k 


Pytes using a 2314 disk pack. This meant that the NPS 


42 





computer could not find the beginning of the second and sub- 
sequent data blocks containing the program. 

DC OuEcetiveraction. eWhen fanally diapnosed, 
fais problem was.corrected by specifying a track overflow 
Meti1on. This option allows data blocksato ‘spill over' onto 
majacent tracks of the 2314 disk pack. Use of this feature 
ms dangerous in that the data overflows onto the adjacent 
track no matter who the adjacent track belongs to or what 
me contains. This requires using a disk pack with enough 
room on it so that the entire program can be assigned 
contiguous tracks. 

Seeernconpatable instructions. 

Peeve wsSant thes track overtlow option, the program 
was successfully copied into temporary direct access storage 
and a load module was created using standard IBM supplied 
Meility programs. At this point the program was in execut- 
able form. The first attempt to use the program was made 
using a bench mark program supplied by ATI which was to be 
used to ensure proper operation. When this program was 
meaaed, the object deck was assembled properly. However, 
Maen the simulation itself was attempted, the NPS computer 
promptly quit and listed an error message which in idiomatic 
translation means 'I don't understand what you are telling 
feo dG and so | quit’. A dump of the core locations con- 
manuing the offending instruction revealed the presence of 
a number of instructions which were implemented for the IBM 


370 and are not executable (understandable to) by an IBM 360. 
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Communication with ATI disclosed that the affected portion 
of the Simulation Program was written in IBM 370 Assembly 
Language and that there were ten of these instructions. 
After further investigation by the author, a similar 360 
instruction was selected that was presumed to operate prop- 
erly for the majority of cases. This new instruction word 
mas placed in the ten locations using a local utility 
program. On completion, the bench mark program was success- 
fully executed and comparison with the version run at ATI 


showed that identical results were obtained. 
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APPENDIX B 


feorrwucli TONS FOR USE OF THE APSS INSTRUCTION LEVEL SIMULATOR 
ON THE IBM 360/67 OPERATING SYSTEM AT THE NAVAL POSTGRADUATE 
SCHOOL 


A. For general and detailed information on how to use the 
IBM 360/67 OS, the reader is referred to the W. R. CHURCH 
COMPUTER CENTER USER'S MANUAL. 

ine AbAG Preerammane Support System Instwhepion Level 
Pmubator currently resides in load module form on Cell 5 
Meme 2o2tevata Cell under file name Sl4S40PBF.APSS. The 
steps recommended for use of this program are: 

iment Meemany cCSiLTed | PLrOgramminmulnG nssenD lye sbanacuace 
mollowing the directions of ATAC PRINCIPLES OF OPERATION 
(VOLUME I). 

2. Include desired simulated Input/Output operations 
following the directions of ATAC APSS INSTRUCTION SIMULATOR 
(VOLUME 2, Section IV). For this purpose it may be neces- 
Sary to replace supplied FORTRAN dummy subroutines with sub- 
routines written to simulate data generation and interrupt 
mening. Ihe Job Control Language sequence required to 
mecomplish this replacement is shown in paragraph C. 

Seine lude desired Simulator Control instructions 
following the directions of ATAC APSS INSTRUCTION SIMULATOR 


mOLUME 2, Section IV). 
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ieerepare the coded version Of lige2, and 3 above fol- 
lowing the directions of ATAC APSS ASSEMBLER/ LOADER 
meuuMie 2, Section II1) for precise _ formats. 

SCONSTIUCE MCG anput cand deck for the desired sjob 
mciluding requiredvcontrol features in accordance with ,ATAC 
meso JOB MANAGEMENT (VOLUME 2, Section I). 

On completion of 2 through 5, the user will have a card 
deck the first card of which is mOpeand the last card is 
SEND. The following IBM JCL cards are required in order to 
SeecuLe the program. 


(STANDARD JOB CARD. SEE USER'S MANUAL) 

TEST EXEC PGM=APPS,REGION=250K 

STEPLIB DD DSN=S1434.PDF.APSS,UNIT=2321, 
DISP=(OLD, KEEP) , VOL=SER=CEL005, DCB=BLKSIZE=400 

FTO6FO01 DD SYSOUT=A, DCB=RECFM=FBA, LRECL=133, 
BLKSIZE=3325S) 

FTOSFOG1 BD DDNAME=SYSTI} 

FTOSF001 DD UNIT=SYSDA, SPACE=(CYL,1) 

FTO9FO01 DD UNIT=SYSDA,SPACE=(CYL, (7,2)), 
DCB= (RECFM=VBS, BLKSIZE=7180, LRECL=92) 

FT1OFO01 DD UNIT=SYSDA,SPACE=(CYL,(7,2)), 
DCB=(RECFM=VBS , BLKSIZE=4 204, LRECL=4 2, BUFNO=1) 

FT20F001 DD UNIT=SYSDA,SPACE=(CYL, (7,2)), 
DCB= (BUFNO=1, RECEM=VBS, BLKSIZE=2004, LRECL=500) 

SYSPRINT DD SYSOUT=A, DCB=(RECFM=FBA, LRECL=133, 
BLKSIZE=3325S) 

Sh4c BY) 


~~ SS SS SS a Se a Oe Oe ee 
~~, ey “TMM, TM, Me, MT, Me, Me, TM, MM 


f/} 
cs 
© 
Ua 


* 


APSS PROGRAM 


$END 
/* 
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B. The source tape for the APSS Instruction level simulator 
Poe ema tieey. kh. CHURCH COMPUTER CENTER under the name 
Zeit tvs”. Filles 1 and 2 on this tape contain instructions 
for the use of the entire APSS system. The instructions 
Pr eine GemnOot nevenchaL appieicable for use and retrieving 
Mime tes iSsneeseonsidered desirable. if retrieval is 
Peeant cd eusce thew ebMeutrlity program TEBGENER. For further 
gnformation refer to the USER'S MANUAL. 

enon eCOle wm meric CD )cen module form of the APSS Instruc- 
Peoneuevelwoimulater. ihe following steps are required to 
prepare it for use in the event of failure to the Data Cell. 
PwemeOt tiese Steps may be accomplished in one JOB if 
desired. The total CPU time required is approximately 10 
seconds. 

ite Ovemtne Ob,ectemodules from the tape to disk 
storage. ENSURE the actual disk pack used has sufficient 
Prablocated room for contiguous assignment of tracks. The 
poyect module file has a block length of 13K bytes and the 
track overflow feature has to be selected in order to be 
Mees co saccess the new disk file successfully. The track 
overflow feature allows data blocks to overflow from one 
track onto the adjacent track, whether it is assigned to 
you or not. The normally used USER disks will not usually 
have sufficient space available. Special arrangements 
Sioulda be Mace, Derore the file 1s moved. The required JCL 


Statements are as follows. 
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// (STANDARD JOB CARD) 
/{ EXEC PGM=IEHMOVE,REGION=150K 
// SYSPRINT DD SYSOUT=A 
// SYSUT1 DD UNIT=SYSDA,SPACE=(TRK, (40), ,CONTIG) 
(1) // DA1 DD UNIT=2314 ,DSN=S1434.ATI.PDF1,SPACE=(TRK, 
ii feceeoelo ye. CONTIG), 
me) // DCB=RECFMST, 
le DISP= (NEW, KEEP) , VCL=SER=PAT003, 
Pi LABEL=RETPD=060 
// TITAPE DD UNIT=(2400, , DEFER) ,DISP=NEW, PASS), 
yh LABEL=(3,SL,,IN), DCB=(DEN=2,BLKSIZE=800, 
Va LRECL=80 , RECFM=FB) , VOL=SER=ATI105, 
yi DSN=ATI.TEMP.APSS1 
// SYSIN DD * 
(3) COPY PDS=ATI.TEMP.APSS1,TO=2314=PAT003,FROMDD=TITAPE, 
TODD=DA1, FROM=2400=(ATI105,3), 
; RENAME=S1434.ATI.PDF1 
ras 


Note mioueialvemianecs selected must be consistent 
throughout. 

Novem tw Intomcelcctcetme track overflow option. 

Neotess: the asterisk (*) in column 72 is 
required as an indication that another 
Siias1OMlOvSavotch 1S agecontinvation 
ot the COPY PDS card. 

2. [he object modules now on the disk must be built 
into a load module and stored in a permanent file. The Data 
Mel eiSeadequate for this purpose. The following JCL cards 
will accomplish this step. 

BUILD EXEC PGM=IEWL,REGION=100K,PARM='OVLY, 

Meee t List! 

SYSPRINT DD SYSOUT=A 
Pere WSN=-S1434.AT].PDFL,UNIT=2314, 

WOh] GER-PAT0N03,DISP=SHR 
Sete wD DSNANE=SYS1.FORTLIB,DISP=SHR 
Se ornov DD Dowame—-51434.PDF.APSS,UNIT=2321, 

Wotan -ceel0S, DISP=(NEW, KEEP), 


LABEL=RETPD=060, 
SPACE=(TRK, (300, 20,2) ,RLSE) 


mA MAL MAL MA MS MRL BL MRL MRL 
™~™, Ma Ta Me, MA, MM, MM, 
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Lil 
[/ 
// 


SYSUT1 DD UNIT=SYSDA,SPACE=(TRK, (19,19),, 
CONTIG) , SEP=SYSLMOD 

SYSLIN DD * 

INCLUDE LIBRARY (PMUDL) 

CHANGE MSIM(IHESAPD) 

INCLUDE LIBRARY (LIBR2) 

INCLUDE LIBRARY (ATACIO) 

INCLUDE LIBRARY (MSIM4,ASEM5,SIMU11,SIMTR1,SIMIO1, GUL) 

INCLUDE LIBRARY (XPLMON) 

OVERLAY Al 

INSERT MSIMUL, *MSIMULA, IHENTRY , IHESAP 

INSERT MINT, IN, OUT 

INSERT IHEDBN, IHEXID 

INSERT IHEBSM, IHECSM 

INSERT Iff[EBSK, IHEIOX, IHEIOP, IHEDID, IHEDOB 

INSERT IHEDIB, IHEDOA, IHEIOB 

INSERT IHEIOA, IHEOCL 

INSERT IHEBSD, IHEBSE 

INSERT IHEJXS 

INSERT IHEOSD, IHEOST, IHEBSI 

INSERT IHEVPF, IHEDMA, IHEVFB 

INSERT IHEDNC, IHEVFD, IHEVFA, IHEVPD, IHEVPB, IHEVSC 

INSERT IHEVSD, IHEVFE, IHEDCN, IHEUPB 

INSERT IHEVFC, IHEVPE, IHEVPG, IHEVQB, IHEVQC 

INSERT IHEABN, IHEIOD, IHEIOF, IHEPRT, IHEVQA, IHESPRT 

INSERT IHEBEG, IHEERR, IHESIZ ° 

INSERT MISET 

OVERLAY Al 

INSERT ASSEM, REWIND, REW72,DSKOUT,CARDIN,DISKIN, 
ERPRT , PRIADD 

INSERT WRDATA, PRICOM, PRINOR,WRITEX, REFTIT,PREF,ERTIT 

OVERLAY Al 

INSERT PARMRD, PRESIM 

OVERLAY Al 

INSERT SMLTR,NRMTRM, STRTSM, TPAGE 

INSERT RDCRD,ABNRMT , ARTHEP , TRACE, HGRAM, HGRAMI , HGRAMS 

INSERT IOINIT,ACTIVE,SIMTIM 

INSERT DEVDTA,ACT,PRGTIM, RAND, TIME, INT, RAND, DEADT, DEBUG 

OVERLAY A2 

INSERT LEVEL, IOR,SPLIT,DMAIOI , DMAIOA, DMAIOD,RIOIO 

INSERT REMACT,DMA,DMATM,RIO,RIOTM,RIOINT, INTOLY , DMAINT 

INSERT DTRAN, PUTACT, RANDOM 

OVERLAY Al 

INSERT HGPRNT 

OVERLAY Al 

INSERT LOADER 

OVERLAY Al 

INSERT LINK, ENTEXT, SLLH 

OVERLAY Al 

INSERT PLATAC, IOPACK 

OVERLAY S$OBJECT (REGION) 


49 





iitomnt OBJECT, INIT, LIB, RCALPH,RCHEX,RCINT 
INSERT MDATE 
OVERLAY $DUMP (REGION) 
INSERT SMDUMP, PAGE 
ENTRY MAIN 
NAME APSS 

/ % 

3. When the load modules have been built onto the Data 
Cell, the disk file is no longer required and should be 
Scravewed. the following JCL, cards will accomplish this. 

// EXEC PGM=IEHPROGM 

// DDI DD DISP=(OLD,DELETE) , UNIT-2314,VOL=SER=PAT003 

WyeeoroPRINT DD SYSOUT=A 

j/ otolNnepe 

SCRATCH DSNAME=S1434.ATI.PDF1 
/ * 

Sess Chto Pegi, sem loodermoatle ison the Data Celt 
Pueethnere are J] memory 20cations which are required to be 
migivilaually Changed tneerder fom the program to execute: 
See APPENDIX A for detailed reasons why. Changing the 
memory locations requires the use of a utility progran 
iV cmMunder close COiLrOleOr tne sCOlplLe: Center. personne! 
for obvious reasons. The locations, current contents and 
meduired new GContents are listed below. For assistance in 
making this change, see the Director of Operations at the 
Sompluce: Center. 


All numbers are hexidecimal 
NAME APSS SMLTR 


Pocatron Old Contents New Contents 
IES8e BDSS 4950 
1166 Biss 4950 
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Location OkeeContents New Contents 
1176 BD53 4950 
iy) BD53 4950 
LIE BD5S3 4950 
1196 BD53 4950 
11A2 BD53 4950 
1200 BD53 4950 
1208 BD53 4950 
IL Le! BD53 4950 
PSBZ ASE B39B 


elie wh ooeinsctructlonebevyel oimulacor 1S mow ready 


HO, use. 


C. A subroutine may be substituted for any of the dummy 
Input/Output subroutines supplied very simply. The first’ 
step is to write and debug the new subroutine. Care must 
be taken to ensure that the new subroutine has the same 
name and identical calling arguments as the old one. Once 
the subroutine is fully ready for use, it may replace the 
original one as shown below. It must be pointed out that 
Piecmuiic sbepers decomplashed, the only way to retrieve the 
Original subroutine is to build the entire system over again. 

// (STANDARD JOB CARD) 

(pe REC EORTC 

(7 eee oo LN, DD = 

(FORTRAN SOURCE DECK) 
/ * 
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Poa iaitieorolLMOD DD DSNAME=S1434.PDF.PASS, DISP=(OLD, KEEP), 
// UNIT =2321,VOL=SER=CELOO05 
Hy eciNKSYSIN DD * 


NAME (SUBROUTINE NAME) (R) 


a2 





APPENDIX C 
PES EM DIGITAL PROCESSOR PROGRAM LISTING 


iicwiurmevrorcal Processor Program (2WPP)" 1s the ample- 
lett aeton on etne algorithn presented invChapter [ilvin ALAC 
Assembly Language. For ease in examination of the coding 
used, a listing of the instructions used in this program 
and their functions is included in Appendix E. 
ioe rocram Listing. Ihe program as listed by the APSS 
Assembler and provides the following information. 

CCOLONI inc firs tecellmnuliSesmeene ins eruction 
Jocation relative to the start of the program in hexadecimal. 
mnie program actual entry address is location 10 so that in 
exces simulator output listaimes, locations appear 10,¢ Ipialyeaiouene 
than they do on the program listing. 

b. Columns 2. The second and third columns contain 
Ene actual hexidecimal Contents of the instruction locations. 
Simele Word instruceons do mot have an entry in the third 
column. 

Gow Une Tinagc sous ticem pit Ppromram Cards 1S reproduced. 
Lines that have a period occuring as the first entry are 
comment cards. 

d. Fhe sequence number of the cards is simply the order 
in which they occurred in the data stream. 

2. The program has been written as one entity in order to 


Pevwemiminiane by SUbrouUtines and paying the associated time 
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overhead. The different functions of the program have been 
separated by blank comment cards for ease of reference. The 
basic function of each instruction has been included in the 
Panmencsmmladc in the margin adjacent to the instruction. 

3. Beginning with the location HADRS, the following cards 
establish data storage requirements and represent the only 
meanSeOnwimterpretimm@eche results ofthe sim@lation. The 
Samulator automatically dumps the entire ATAC core storage 
Mmeonsheechon OL theesimulation., By referencing the loga- 
tions of the various variables andudata, thewcontents at 
ene completion of the simulation run may be determined. A 
summary of the meaning of the data storage locations is 
maven below. 

Variable lnipenpre Gaicom 

HADRS Used tomepecifyelana¢smof the instruction 
versus time Histogram 

TRACK tee Otte pulSe train file. Individual: 
files begin every 10,¢ FOCALIONS = IN tenpue— 
pane ee. re ie noted values is contained in 
Seccrone the Cel a7 5., 

ER Refequency inpuc tile. ihe contents of 
EFILE locations show the mapping conducted 
from frequency information to the proper 
pulse strains. 

Eas The €contace filesentry file. Locations 
Comntespona fomtic puls@ train file and are 


in the same order. The contents of even 
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Variable interpreta tation 
words of the file contain the entry location 
GEeltnemeonrespondineapulse train tile. 
OFLO The overflow storage area. 
OFLOW Variable used to point to the next available 
locatron ei nerhicwover Slow tated. 
NEXTF Variable Used tow point stomthe last@pulse 
train sre iemlvateds 
Pee bine Used to defines the extent of TFILE fou 
the purpeses wt inetializing the TFILE 
TOcations to the proper value by reading in 
data cards. 
DATA LOCA ttomemtsSead  tOmseore Inputepmlse data. 
DATAR Usicdetoadehinepthemextent of DATA for tiie 
purposes Of initializing the DATA locations 
to the data values by reading in data cards. 
4. The simulator output. 
ihe pt Tit eOor stheCwContenusuormeadcen Of the Lo ceneral 
registers, the condition code register, the input/output 
-rcurupt Masks and the Simulated computer time elapsed 
Since the beginning of the simulation can be provided. This 
print out must be specifically requested with a Trace 
PeCducict we leo 1S etunmMead son and Off within the program to 
avotaeda line Of Output for each and every instruction 
SS oateedshy etic Simulated computer. The feature 15 very 


useful for checking the program flow. 


Se 








b. The simulation histogram of time expended by the 
Mmistructions at G€ach real location (10 + the program listing 
location). This feature shows which instructions take more 
fime in execution than others. By noting which instructions 
wre instruction Sequences take abnormal amounts of time, the 
mrogram flow can be changed to improve the overall effici- 
ency. Additionally this output shows those instructions 
which should be looked at with the possibility of redefining 
a new microcode instruction or modifying the hardware con- 
figuration in order to improve the time required for 
execution of the program. 

ew COre duUNpe mw lseoueemt 1s the primary means of 
verifying proper operation of the program. Variable 
MecacionsS in the program listing are given in actual core 
locations if the variable is in Random Access Memory (RAM) 


or Read Mostly Memory (RMM). 
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APPENDIX D 
THE JFM DIGITAL PROCESSOR PROGRAM DETAILED DESCRIPTION 


A. PROGRAMMING 

The coding used for the ATAC Assembly Language Program 
was developed by encoding each block of the algorithm flow 
chart of Figure 1 separately and then linking the blocks 
mosciner.  lhese coding bllocks are described below. 

1. LOAD DATA FROM DATA STORAGE. The data loading 
Seetation 1S Critical in that it requires considerable 
-ceue@on timc and 1S reqismeamror everyepulse. The most 
Peterlee load Instruction 1s ene REIURN@ which is designed 
to return stacked values. This instruction automatically 
fetches a sequence of words and decrements the address 
foitieera il Using this =tistruetion in this application, 
prior and post operations to the address pointer are 
Paoircdsinisorder to kKeepeche poinwer at the proper location. 
Micwscaquence required 15 as follows. The actual instructions 
ane indicated. 

TOP ADD IS,0,4 Add 4 to the address pointer in 

regiceer 0. 
“ADD LS ase! Merete emerogram Control Pulse Counter 
Bice 2,DCNE Go stop simulation if enough pulses 


have becmeane lyzed. 
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hel DR,O,4 Starting at the address in register 0 
decrement the address and fetch the 
contents of the address. Repeat until 
4 values have been fetched. The re- 
turned values are stored in registers 
1-4 in the order they appear in memory. 
On completion of the operation, the 
address pointer value is 4 less than 
eee ne DECOR nian o 

eet TVA ene TO test memescCONntenteseor register 1 

an operation has to be performed to set the, condition code 
meeister. lhe operation selected 15 to move the Contents 
imc tIster 5 fOr futune Operations. Kegyster 1 will later 
be used as a general zero contentS and indéxing register. 

POR “kh, 5, 1 Load the contents of register 1 into 
rerester o. 

NCOs . LOr If the contents of register 5 (frequency) 
are Zero, Uranen to the peginnine (IOP) 
and access the same data block again. 
This sets up the no input data idle 
NGO. 

3. ZERO DATA LOCATIONS, INCREMENT ADDRESS POINTER. The 

frequency word location in the data memory is set to zero 


Mdleatine stwewaata Lor that pulse Nas been read. The address 
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mpotiter in register 0 1s incremented to point at the location 
of the data words of the next pulse. 

DR IS jae, 0 Load 0 into register 1 

SiR RX,1,0,1 Store the Mant ents of register 1 in the 
location whose address is the sum of 
registers 0 and l. 

Abe 1S 2054 Add 4 to the contents of register 0. 
ce teriemacttemadLe pointer for the Wext 
pulses 

ae) PEND CONTACT “Pive NuMpER = Iie wt requency worden 
MmoloterT 3 1S USed aSmae relative address for entry into the 
frequency file. The frequency word is copied into register 
7 SO that it can be saved imethe event multiple files exist 
mom this freduency. 

Pee RR. 705 LOddwinewGontentsme: register S/amco 

TeCaisier aye 

mmo R DxW6,EFILE,7/— Load=reciaster 6 wath the contents of 
the location whose address is the sum 
of EFILE, the base address of the fre- 
quene, files and resister 7. 

Teo tbcarmeholss If oO tile 2e6nrmgenss frequency has 
Deciecotablished, #he contents of the £requency file is 0. 
Peimeane lm On COnditlion zero inStruction 1S used. 

BReOn aL, OfLARL ime Qunenley a wOraesrsezero, branch to 

“STARI where a file will be created. 
fee!) CONPACGCEAEILE ENGRY ADDRESS. This anstruction 


1s executed only if the frequency file already exists. The 
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Povenre byte ot the frequency fille entry is the relative 
Pfaress of the proper contact file entry. This file Hee 
two words, the second word is the pulse train file entry 
address. The coding to achieve the linking is shown below. 
At the start of the sequence, the frequency file contents 
are in register 6. 
LDLB R,14,6 Load the lower byte of register 6 
mneo Tegmwereim.: 14 
LDR DX,8,TFILE,14 Load register 8 With the contents 
of the address formed by adding the 
contents of register 14 and the value 
Ope Jes 
7. PW = FILE PW? The most "efficient coding of Chis 
algorithm step depends on the relative frequency that the 
momdition is met. If tht’ avera@ge condition is that the 
incoming pulse width at this point will be equal to the 
file pulse width, the most efficient coding will return the 
zirst two words of the pulse train file together. If the 
branch is to be taken, the most efficient coding will 
return just the single pulse width word. The latter con- 
dition was chosen. Using the entry address contained in 
mesgister §, the first pulse width of the file is returned, 
Sompared to the incoming pulse Width and the appropriate 
Seranch taken. The ATAC instruction set contains a compare 
between limits instruction which would perform the desired 
operation but it does not have a mode that is compatible 


With the program. The use of the ATAC instruction requires 
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enough manipulation of the contents of registers so that 


Pemmequires more execution time than directly coding the 


individual steps. 
follows. 


ihe RX. ode 


Pee R , 1S y.9 


SUB R,13,4 


BRCS P,TEST 


Pecmen, 1s, 15 


Geol Serres, 13,4 


BRCS P, NTRY 


The coding for this operation is as 


Load register 9 with the contents of 
the address formed by adding register 
8 and register ] (register 1 still is 90) 
Load register 13 with the contents of 
register 9, This is done to eliminate 
the requirements to load register 9 
again. 

Subtract the contents of register 4 
(incoming pulse width) from the con- 
temic et reels tcmel>s ands placemenc 
results in register 13. 

If the result is positive branch to 
the positive test (skip the following 
1S SR Ueehon ue 

Load the arithmetic complement of 
YEG@Ister loeeivwommocceer 15. 

Compare the magnitude of the pulse 
width difference in register 13 with 
the threshhold value 4. 

Tf ENE wCONLeH tse Ormmee ister 15 15 
Creare weenie ecm bmaAlIch 1S take 
to llocabioneta ies wnaeh 1S the address 
of the beginning of the LOOK FOR NEXT 


FILE routine. 
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Se tino FRI BEEN CALCULATED? The contents of the loca- 
tion following the pulse width word in the pulse train file 
asthe PRI control word. If its contents are positive, a 
PRI has been calculated for the pulse train and the contents 
Sautme word 1s the value of the PRI. JI£ the contents of 
fmemwerdad 1S 0, the pulse train has an unstable PRI. In 
bOtn Of these cases, the program should proceed to update 
micepueesce tragn file data. If the contents of the PRI word 
Wecmless tian Zero, the pulse trawmerile is still building 
Momence precman Should preceed te tne=puinda tile reueine . 

SUE Powe. 1 AddViteomencweontents Of register 1: 
Register 1 will be used an an index 
to step to the next word in the pulse 
tPaaiieeed Ler 

LOR ene, LO, 1,8 Load register 10 with the contents of 
thesaedycssmtormed (by adding the con- 
tenes Of Fecrereroel and 6. 

BRCS N,BUILD Branch von tene became or the build 
fi ler outline Lee tlemecnrLents Of 
regirocer mare mmlecative. 

eel wesw oI NCE LAS! ibs relr sche ine UPDATE route 
now calculates data to be used to update the pulse train 
Ci el ttromeirst Operation 1,eaquires the returm of the next 
ivOewOrdS Of the pulse train file which contain the tune 
fat the last update occurred. The difference between the 
last update time and the current pulse time of arrival con- 


tained in registers 3 and 4 is used to determine whether the 
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PRI will be updated or SPR calculated. The test used is 


whether the upper half of the double precision difference 


of the two times of 


diem | GeezZere. las results in the 


algorithm considering all time intervals greater than 0.8 


seconds as SPR and all intervals less than 0.8 seconds as 


PRI. 

UPDATE ADD IS,1,1 
ER Re eles 
ADI) 1S IST 
TDR RG pa Se 
LDR R,6,2 
EDR hs 


POURS ko, 11 


Scep the indexing register to fetch 
INE Tek WWweneal 

Load the lower half of the time word 
into register 11 

Step the index register agagm 

Load the upper half of the time word 
into register 12 

Load the contents of register 2 into 
register 6 

Load the contents of register 3 into 
register 7. These two instructions 
copy the double precision new time 
into registers 6 and 7 so that the 
aculalwtame Value 1S retained an 
Beo US tCiceszeand © S 

Double precision subtract the contents 
Otwmecisters llsand 12 from the contents 
of registers 6 and 7 and place the 


Pesiltseonerceiusters 6,and 7. 
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icons 1 OP liecnemte sult 1S negatage, diasregard 
this pulse by going to the beginning 
tom cete tlic Next pulse 

BRCS NZ,SPR If the upper half of the double 
DEC TON rTerence aS not zero, 
DaaicimeOneNemoiwemcalculatizon routine. 

LDUB R,13,6 Load the upper byte of register 6 into 
Ee CiSitren 615% 

GP Wee 15,16 Geipareetnewcontents Oo: repaster 13 
with the magnitude 16. If the con- 
tents of reguster 13 is greater than 
ioe the tame iamterval between pulses 
1S greater m@ehane 52.4 milliseconds. 

RCo Ps SPR If the time interval is greater than 
SZ Osi Mitccconds.coecalculate SPR. 

10. UPDATE PRISPWOCAND TIMES .PRivand PW.are updated 

Ppecceecrmining themmwetaiimivelecicemor the new value and the 
Opceee li@ene new valulemis=greater than the old, 1 is added 
COm@re Orarand Vicemmersaamelt therewrs no difference, no 
Sewiecemtceimlade. © [Hemme IS5SeSInmly sioaded 1nto registers 
11 and 12. On completion of this step, the four updated 
pulse train file words are contained in registers 9 through 
dae 

Sup ko, 10 Subtract the contents of register 10 
from the contents of register 6 and 


DIA eCeeGliGwiEesina: iii register .6. 
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BRCS P,INCR 


BRS N,DiEGR 


SORE LDR R,11,2 


LDR R,12,3 


SUB R,4,9 


BRCS P,INCP 


BRCS N,DECP 


If the difference is greater than 
zero, the new time difference is 

ia Geet ee manchsto the INCRease, routine 
which will add one to the old PRI. 

If the difference is less than zero, 
the new time difference is smaller. 
Branch to the DECRease routine which 
Vit olbenaGeteole arom the old PRI. 
Voadetncemem time located in register 
Z into register Ill 

bOaGmtiomuope Galiaiteotethesnew time 
LOCA dmiiGeCre tone seintonneoaster 12. 
lien GE Memeomuen iS wot oreoas tei. 9 
from the contents of register 4 and 
plLAaGemencwue SUT omInirTepister 42qe This 
subtracts the old pulse width from the 
new pulse width 

If the result iS positive, the new 
pulse width iS ereater. Branch to the 
increase pulse width routine (INCP) 
ech Wilileaddistomtine Old pulse 
width. 

If the result 1S negative, the new 
pulse width is smaller. Branch to the 
decrease pulse width routine (DECP) 
Wiltewewilieenbtract | from the old 


pulse width. 
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tis STORE UPDATED DATA. The four updated file words 
Sree placed in Wemory uSime@wa stack instruction. Aiter the 
stack operation, the program branches unconditionally to the 
wecinnaing LO cet the nexe pulse 
eeck  SIJK DR, 8,4 Store the next four registers (9-12) 
into memory locations beginning with 
the address contained in RegEster 16, 
Add=45°t6 the coments Ofmmeo1 Sterns 
Ol Conn uer Lon. 
Pies U, LO? Unconditionally branch to the beginning 
(TOR) Stemeete these seamusise. 
iZ. -STARYT PILE. @ihemstanteeemomeme “(Sikes 
entered when the frequency file word is zero indicating that 
MOwoulSe trail file e6xXasts™ cor tic fa@cucncy detected. It 
is also entered after linking information has been prepared 
moeme the first file waS@WeC the correct one. The routine 
Sepencs Lor the Next Empty pilse train c£ile gwd then enters 
‘linking information into the frequency and contact files. 
ipewddata 1S then loadediamto the firstGemtery Laeation with 
—lomstored in the PREMwond) indaeat inceencanuilse mt oats 
MileceHas been receivedueeihe loeeronmhieireaseused to 
point to the last file activated so that the files are 
activated in sequence. A mask is used so that the values 
in NEXTF will always be modulus 64. Files are cleared by 
the main classification computer. 


Seem LOR fl 13,-193 Load m¥edulo 64 mask into register 13 
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LDR D,6,NEXTF 


ANOT ADD IS,6,2 


AND R,6,13 


MR De, 8. TFILE. 6 


BRCS NZ,ANOT 


STR D,6,NEXTF 


STR DX,5,TFILE, 6 


boaew time contents of location NEXT? 
Mieor recs ter 6 
Register 6 contains the last file 
used. Add 2 to get the address of 
Ene Next wile entry . 
Form the logical product of the contents 
of registers 6 and 13 and place the 
results in register 0. = le = ecomuecntus 
Of Meciste? oe dre no Modulomoae 
Load in register 8 the contents of 
the location whose dddvess ascethe sum 
Gt ETRE 2n¢mreorstei oe wiliior Th omeeie 
address of the frequency associated 
WLU Che TicKkc Pulse train 
If the frequency value 15 not zero, 
branch to location ANOT which is the 
MiStruceLoneco Increment Neem. 
Store the contents of register 6 in 
location NEXTF. This stores the 
Pomicen Valen loOnmene Pulse train tile 
which 1s now being Setav awed. 
StOmemene CONbenes of recister 5S in 
the location whose address is the sum 
Of TERTLE and@regserer 6. This speci- 
fies the frequency of the pulse train 
file and also indicates that the file 


LS DUSY - 
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ADD 


Lo cel Add one to the contents of register 6. 


Pree DX,8,1TFILE,6 skhoad register 8 with the contents of 


STR 


LDR 


LDR 


LDR 


LDR 


STK 


the location whose address is the sum 
of TFILE and register 6. This is 
the second word of the contact file 
entry and is the actual address of the 
first word of the pulse train file. 

DX, GO, EF iligy Store the contents of repismer 6 in 
the location whose address is the sum 
Om rbiLE and recister /..0 has sets the 
link between the frequency word and 


the pulse train file entry address. 


R,9,4 Load register 9 with the contents of 
repistereai(pudseevid thos 
FSee = 16 Load register 10 with the number -16 


(pulse counter) 

Reel Load register 11 with the contents of 
register 2 (time, lower half) 

Repel Zac 5 Load register 12 with the contents of 
FecistersomutIme upper halt jo the 
pulse data is contained in registers 
Om eh otic lhaeleZe 

DR,8,4 Store the contents of registers 9 
through 12 in memory locations beginning 


Wier tiem aindtesomenmetanoea ean regastcr &. 


mies U,TOP Branch unconditionally to the beginning 


(TOP) to get the next pulse. 
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13. WRONG FILE. The routine NTRY is entered when the 
pulse width of the previous pulse train file was not within 
Timits. This routine establishes linking to the next pulse 
train file by taking the upper half of the frequency file 
word and placing it into register 7 and branching to FLINK 
if a file exists or branching to routine STAR2 if another 
file has not been started. FLINK is the execution step 
mideecntecrs the frequency tile. The contents of repister 7 
Meow be tne Overi low location in the file. STARZ esta- 
blishes the initial linking for establishing multiple files 
Pomeeme Same Lrequency. 

Pere LDUB R,14,6 Load register 14 with the upper byte 
of the contents of register 6. | 

BeGoee, of ARZ Peernewunpes byte oF recister 6 was 

Zero, anotier pulse train file has 
not been established for this frequency. 
PranvenetO olAneze and establish one. 

LDR R,7,14 Load register 7 with the contents of 

bees ter =lad’ 

BRCS U.FLINK Unconditionally branch to FLINK. FLINK 

is the llocation of the FIND CONTACT 
FILE NR sequence. 

14. START FILE (2) The routine STAR2 is entered when 
two Or more files are eStablished for one frequency. On 
TiMwemmee slower byte or repister 6 always contains the 
relative address of the last frequency file accessed and 


Miemuppenr byte 1S always zero indicating that there are no 
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Tieemiieres ii tie Ciajin. STARZ sets the upper half of this 
mem word equal to the address of the next entry for this 
frequency and stores it in the word's proper location. 
STARZ then branches to routine STAR1 which finds the next 
empty pulse train file. 

LDR D,14,0FLOW Load register 14 with the contents of 
leocation@OFrley. “OFLOW contains the 
relative address of the next available 
frequency file overflow location. 

ROB R,13,14 Load register 13 with the contents of 
reguseer Adit Che tnesby bes wemenanced:, 
This places the overflow location 
relative address in the upper byte of 
register 14, 

POD O11. Add the contents of register 13 to the 
contents of register 6 and place the 
TresultSein Terpister 6. this assembles 
a word im register 6 which consists of 
two separate bytes. the lower byte 
1s the relative address of the frequency 
file entry that lead to the incorrect 
pulse tiyain fake. Hiesuppem byte is 
the relative address of the frequency 
fi lementryaivnerehs wiidieakead to the 
COnmrect epulsemenaime ft) leeiorsthemulse 


width ofethemecurment pulse. 
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STR eoweiees,/ orore the contents of register 6 
MiechestOcattron whose address 15 the 
SUM Ouse | beeand repaister 7. 

EDR R,7,14 Put puower relative address into 
Yeors tee, 

noe to, le, 1 MOGeOne = tO Cle “Contents Of register 14 

ew rreor GOW Store the contents of register 14 in 

location OFLOW. The contents of OFLOW 
(ove VOMMmeo Lie Nextwavarlable @ver - 
flow location. 

bee U;STARI Unconditionally branch to STARI1 where 

tine =a oemierain £1le will be 
2G tigated. 

Pere Ouard PULSES *TOVCALCUBRIE PRI? Routine BUILD is 
Smecrecamyien it has been determined that the PRI word of a 
femeceetrain tile 1s léss than” zero indicating the file is 
Bapeleouaidaing. The number one 15 added to the negative 
fPomisemeouUnter and it 158 tested. If the counter is now zero, 
Meemecurrent pulse is the last one needed to calculate the 
Momdicgetie program branches to the calculate PRI routine 
feet the Counter 15 still népative, the program then 
Continues to add the pulse data to the pulse train file. 

mao eet, 10,1 idm tOm tuner Coneents Of repister 10. 

breoncl, RI If the counter is now zero, branch 

tomenewealculates snr 1routine= (PRI). 

Pee De NFORMATION TO THE DATA SET. The pulse counter 


is stored. The proper location of the current pulse data 
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is determined by adding 16 to the pulse counter value. The 


result is the number of the pulse in the current pulse 


fee 11). 


This number is then multiplied by 4 (shifted left 


2 bits) to determine the address of the first data word 


relative to the start of the file. The data is then stored 


azter checking to ensure it is valid data. 


STR RX,10,8,1 


ADD 


SHS 


ADD 


LDR 


IS,10,16 


LL,10,2 


R,8,10 


Store PheMecontentseoterecaster 10 in 
the location whose address is the sum 
of registers 8 and 1. Register & 
contams the address of the first loca- 
t10n OL the pulsecmtrain £lle. Register 
1 contains the number 1. The pulse 
counter is stored in the second word 
Gpethesti le, 

Add 16 to the contents of register 10. 
ie Gesiht 1S thesnunber got the slot 
LOr thes cur renmiuaypmhsem@data . 

Site the Conmeeiec on soouster WM Lert 
SWORD TLCS eI ilsminuleci plies the contents 
of register 10 by 4. 

Add the contents of register 10 to the 
contents of register 8 and place the 
results in r@@ister 8. This forms «the 
actual address of the first word of the 
l6Ocat tons ehoOmebiemGcunmrentspulse . 

Load register 5 with the contents of 


MEG LStCl oO. 





ADD 158,5,-2 sumeract 2 from the contents of 
Fegister-o. This tommemthe addwess 
of the time of arrival of the previous 
PiRESeROLNthewecuTrrent pulsestrain: 

LDRM DR,5,2 Load registers 6 and 7 with the contents 
of two adjacent locations whose addresses 
Sear t demene. CONtente 0.) ogi cher 25. 

BOR Rs bh 2 Load register 11 with the contents of 
register 2. 

hone IR ,LZ. 5 Load register 12 with the contents of 
Fepwuster ole Ine  Ommiaaer ty alelO1: 
the’ Cupgent. pul wna bccn COpmeds Into 
registers 11 and 12. 

WOW De 5125 0 Double precision subtract the contents 
oferecnisters beandey from the contents 
OrerepiSterse2 elaeos and’ place the 
results in registers 2 and 3. This 
puts the time intervads betweem pulses 
i ReGrstews ec samc. 

Bes NN, 1 OP If the time difference is negative, 
branch to the beginning disregarding 
Che Cupmemeulse. 

BReo. Z,GOOD LibthesuppereialGwotesthestime difference 
1s) zero, the time 2s considered a PRI 
anc heaie 1S ERNeGmTOn 1S enot 


executed. 
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MOR 1S, 2,0 Load register 2 with the number 0. 
This will cause zero to be stored for 
the interval between pulses. 

GOOD Le Komal). 2 Load the single precision value of 
the interval between pulses into 
fegister O0=rrom register 2. 

LDR R,9,4 Load register 9 with the contents of 
register 4 (pulse width). The contents 
Of Teseisters 9-12 15 now the data 
associated with the current pulse. 

Shen. 6 4 Store the contents of registers 9-12 
in locations whose addresses start with 
the contents of register 8. Add 4 to 
the contents of register &. 

Reo. U, TOP Unconditionally branch to the beginning 
(TOP F=t0 svete ncmmlc Kt pulse. 

17. CALCULATE AVERAGE PRI AND PW. When the seventeenth 

pulse of a pulse train file is received the program enters 
PRI whien calculates the average PRI and pulse width of the 
Mulse train file. The time anterval between the last two 
pulses is calculated and then all sixteen intervals and the 
last sixteen pulse widths are added. The program then 
branches to PEND where the sums are divided by 16 to form 
the average. 

imeiok oe is,1,-15 Load Tegisterelewtrimtneanumber 15. 


Mrrs 1S Uscaumas eameoun ter 
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ea 


ADD 


REE 


LDR 


LDR 


IS,8,64 


PeUtoeek .O 11 


bates N, LOP 


LDR 


IS,11,0 


POR, 12,6 


ADD R,4,9 


Add 64 to the contents of register 8. 
Register 6 MOWwecontains the required 
address for the RETURN instruction to 
return the last pulse train file data 
Sete 

Load registers 9-12 with the last data 
set. 

Foad@recrster"6 with the contents of 
register 2. 

Load register 7 with the contents of 
Resmster Soawelne tame Of arrival of 
ChescimneieSeeES —COp1ed 1nto regis - 
ters, 6 and 7. 

Double precision subtract the contents 
of registers 11 and 12 from the Pontents 
OimecustcrS Gmcmida/. The lasitmtime 
interval is contained in registers 6 
and 7. 

if tie Lame Gditterence is negative, 
disregard this pulse. 

Loads0 into register 11. 

Set ee iestandard. 

Add the contents of register 9 to the 
Contents ot register 4 and place the 
Sum in register 4. This accumulates 


the sum of the pulse widths. 
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DADD R,6,10 


ADD US el sal 


BROS Z., PEND 


ADT Syaae, - 2 


RET DR,8,2 


Ss Oa 


Gor fF , LPS 


CMP 1S,10,-4 


BRCS N,NG 


ABD R,10,12 


Bicol, PZ 


LP3 CMP IS,10,4 


Dowwle precision add tie contents of 
heguoeciceH0eand sllato the contents 

OE MECC OMS tens moOmanGd “7 

Add) 1 to the contents of register l. 
THUS™Steps the counter . 

Prwenmewcounter 1S zero, all pulse data 
has been summed. Branch to the divide 
routine (PEND). 

Setwaddmesis for thiemmext data wet. 
Return the pulse width and the time 
interval from the previous data set. 
Find the difference between this time 
interval and the PRI standard. 

limite Mat tenencewis positive, test 
themresults ait Geer 

Compare the difference,-when negative, 
With the meéemitude -4. 

[mechne “ditterence 1s less than -4,; 
gemthrow thempulse out. 

Restore the Game seiterval value in 
Feguster dtr 

UneGonGataonaldysbranchieto, LP2 te sum 
theedatax 

If the difference between the time 
Miter val Maiomeiewniestandard 1S posi- 
MiMemmeonvarewunecrati ference to the 


magnitude 4. 
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BRCS P,NG 


ADD R,10,12 


BRCS U,LP2 


NG LDR R,10,12 


BRCS U,LP2 


PEND SHD RA,6,4 


SHS RA,4,4 
RET DR,8,2 
AM 18,8, 2 


LDR 1S,11,0 
wR 15,12,0 
STK,DR,8,4 


Preven itdemence 1S greater than» 4, 
£O thmewethe pulse out. 

Restore the time interval value in 
Tecisiter el). 

Uncongeeedemally braneh to LP2 to sum 
the data. 

LopmceemtnomrEanie Interval wath the 
value of the PRI standard. 
Unconditionally branch to LP2 and sum 
the data. 

Doublouprecisaon Shirt richt regis - 
ters 6 and 7 4 bits (divide by 16). 
Shit t Reise enh bts tO 
divide wiyliGe 

ROtupmetiomierrst €lnemot arrival in 
te eae 

Set the address for the second data 
block, 

Load 0 into re@ester 11. 

LoddmUmiitto Tem ester sic . 

Store the contents of registers 9-12 
in locations whose addresses start 
with the contents of register 8. This 
Seores theetimes: required to calculate 
SPR and sets the initial value of 


SER com un 
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Perk Ro, 4 Load register 9 with the contents of 
rectecen 4. 

LDR R,10,6 Load register 10 with the contents of 
register 6. 

EDR R,11,2 Load register 11 with the contents of 
register 2. 

eR R5i253 Load register 12 with the contents of 
register 3. the entire data set is 
Poadeqd wine Tregasters J-12. 

ADD I1S,8,-8 Set the address of the first data block 

of this pulse train. 
STK DR,8,4 Store the contents of registers 9-12 
in locations whose addresses start with 
the contents of register 8. Add 4 to 
the contents of register 8. 

Beak UU, TOP Branch to the beginning to get the next 
OwISe- 

Poop oALCULATE SPR AND UPDATE TIME. Routine SPR is 
entered when a large time difference is noted between the 
current pulse and the previous pulse of a closed pulse 
train file. The time that the last sequential group of 
pulses arrived is accessed and subtracted from the current 
time to determine the antenna scan time. The current time 
aaa@mene oPR are stored in the file. 

Seow 1t5,1,1 Increment address pointer. 
Were kxX,11,8,1 Load repigu@mer! with the contents of 


thewlocativom whese address is the sum 
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of register 8 and register 1. This is 
thesstare of the time that the last 
antenna scan was received. 

AEP IS, 11 Load second half of this time word 


DR RX,12,3,1 HiteOr he Ciesiter lz. 


ORR, O52 Move the new time of arrival into 
PRS eR 7,5 WOT Ile mols GCIs 
DOU —=hur6.,71 Double precision subtract the contents 


Of registers Ll ana WZ trom the contents 


of registers 6 and 7. This is the SPR 


value. 
EUR obs, 154 Set the address pointer 
SiR RAs 2,1,6 Store the new time of arrival 


Ply io. 1.1 

elk  RxX,5,1,8 

oy eS 51, Ll 

Stk) kx,6,1,6 Store the value of the SPR 

ADU TIS, 1,1 

puke RA, 751,65 

pico U,o1LORE Unconditional soranch to STORE and 


LINISh Updating sthe pulse antormation. 
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APPENDIX E 


DATA SIMULATION PROGRAM GLOSSARY OF VARIABLES AND 


FORTRAN LISTING 


DATA SIMULATION PROGRAM GLOSSARY OF VARIABLES 


ANG(I,J) 


ANT (I,J0 


BIN 
BEAM 
BNG (I) 
BNGI 
BW(I,J) 


CON1 


COMNZ 


SEL 
CoBo 
J3u 
D2 


True pointing angle of the antenna on the Jth 
emitter on the Ith platform 
SMeeiieSCaMmGdie OL stile mtn emitter on the 
Ith platform 

EW Sehequenc yy same Suze 

Input value of antenna beam width 

Irue bearing of tine Eth platform 

Computation value of BNG(T) 

Antenna beam attenuation factor of the Jth 
Cie ver re On sUhc timp matron 
COnvVersSTOnyiaectopeused win Converting from 
antenna beam width in degrees to the 

antenna attenuation factor 

COnNVeETrsion factor uUsea sin converting from 
Cransiicter peakepovweneco transmitter total 
effective radiated power 

Input value of the course of the Ith platform 
Input value of the course of the IFM 

Distance to CPA 

Distance to the maximum range at which an 


emitter can be detected 





DEGR 


DEL 


it 1 sed) 
FMAX 
FMIN 
FREQ 


EOCALE 


GAIN 
et 


IF(I,J) 


INEXT (L) 


IPRI(I,J) 


Pe Gee 


Conversion factor to convert angles in 
degrees to angles in radians 

Ditecaierctmine Start of the samulation that 
aiecimeteGomtirst pulse is transmitted or 
that an emitter is within its maximum de- 
téectiGne range 

Transmitted frequency of the Jth emitter on 
Che wien laerorm 

Mos smuneucccivemareaquency ot the ILFM 
Minimum receive frequency of the IFM 

Dipti luemOt ste Operating trequency of 


an emitter 


Conversion factor for converting from Mega 


units to units 

Input value of antenna gain 

COUM@CEN LON MUNDemmomepulses received by 
the IFM 

Selceto, 20) the desired type of frequency 
variation 

Platform subscript of the second or more 
pulses when more than one pulse is received 
Simviseancous ly spy ethos LEM 

Sc lieeeoO) sonmetiecug@coimedecy ne of. PRI 
Variation 

Sc le eeen Mier cClicu@ccimeqatype of output 


power variation 


94 





‘ 


IPW(I,J) 


ISCAN (I,J) 


ITEMP 


JNEXT (L) 


cae MP 


NBIN 


NGAUS 


NMTR(I) 
NOUT (K) 


NPW 


NTIME 


NTIME1 


NTIME2 


Selector for the desired type of pulse width 
Variation | 
Pomc TOnetnemacscired type OLeSecane nace 
Variation 

Platform subscript of the clyrent pulse 
Received DY stne un 

Emitter subscript used for multiply received 
pulses 

EMDtter sSuUbSCYioe of the  enreent pulse 
received by the 1FM 

Bin number of the frequency of the received 
puese 

Number of frequency bins in the bandwidth 

Of tne 
Integer calling argument for the Gaussian 
Distribution Random Number Generator 
Number of emitters on the Ith platform 
Array used to prepare output data for the 
ATAC Simulator 

Integer representation of pulse width as 
required by the ATAC Simulator 

Integer Veprescmuatmoenean SIMULATION Gime 
Power Nhalt Ot  themdouule precision. time 
value required by the ATAC Simulator 
Upper halt On (tiesdoubmesprecision time 


Value TeCquiredp, —tne ,1AG Simulator 
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NUNI 


PI 
PI2 
POWER(I,J) 


POWIJ 
PRI (1,J) 
PRIN 
PSCALE 


RCVBNG 
RCVR 
RECOIG 
Greve 


RDEG 


RMAX 


SCAN 


imeecer calling argument for the Uniform 
Distribution Random Number Generator 
3.1415926 

Pay 2 

lOGCALVerreetive -GrallemlEted power Of “the men 
emitter on the Ith platiorm 

Caren latrOn yalierom POWER Ls. ) 

PRI of the Jth emitter on the Ith platform 
Input value of PRI 

Conversion etactor £Or CONverting trom mle 
Lis LON Lts 

Input value of transmitter power 

SilenlAGion wallne of PW(I,J) 

Pulse width Of thet eemitte, One iic as 
platform 

Range of the Ith platform 

Bearing rromewnien. tne current muse il lames 
Prequeney On suiemiceel Ved pinuce 

Sienal level of fenewreceived pulse 

Direction of Relative Motion 
Commersiqnmiacton onmeonverting sane1co mn 
radians to angles in degrees 

Maximum range at which an emitter can be 
detected 

Subroutine that ealenaitcs an antenna true 
POLNtaAng ang le Vande in tired rection 


Gt encoskinm 
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SCNMIN 


SEND(I,J) 


SIGFRQ(I,J) 
SPDI 

SPDO 

TF 

THETAI 
THRESH 


TINT 


TNEXT(I,J) 


TO 


TOF 


BO CALE 


TSKTS 


TT 


TWOPI 


UPDATE 


VL 


The minimum value of antenna gain for which 
an emitter can be detected by the IFM 

Time of transmission of a pulse from the 
Uti cence l Onewemnce lth plakitorm 

Frequency of the transmitted pulse 

Input speed of the Ith eplattormm 

Input speed of the IFM 

Simulation stop name 

Bearing of the CPA of the Ith platform 
Threshhold sensitivity of the JFM 

Minimum time measurement capability of the 
Sysvem Util tino sincere, 

Time of arrival at the IFM of the next 
pulse transmitted DY “thee (a) CM tte rroreeas 
ith pletion 

SIMA tones teart tine 

Time required for the propogation of a 
pulse from an emitter to the IFM 
Conversion factor for converting from Kilo 
LEME ins 6 Oy bias 

Conversion factor etorecomvertine from knots 
to yards per second 

Calctlation valuc ofeINiE iil.) 

ZENE 

Subroutine that calculates the time of 
arrival of the next pulse from an emitter 


WeMLOCILY Of b1iene 
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VTI 
VX(T) 


VYO 
WSCALE 


a) 
XI 
XO 
Y (I) 
YDS 


YI 
YO 
el) 
ZI 


ZO 


Motal reWative velt@eity of the Ith platform 
Pelaurwemvelocity Of the 1th platform in 
the X direction 

Calculation value of VX(I) 

Velocity of the IFM in the X direction 
Relative velocity ofsthe Ith plattorm in the 
ema 2 Cres eakon) 

Calculation value of VY(1) 

Velocity wOnmLiew lL Malin tine w\ ecdine ct som 
Conversion factor for converting from micro 
units to units 

Relative X position of the Ith platform 
Input eesitionseL Ehe th plactom 

Input Moves e101 of the lrM 

Relative "posttrolmoretnic [ih platiomm 
Conversion factor sronr Converting snocs £0 
yards per second 

Inpug.Y positionvepmenc lth plattormn 

Input Y position geieene IFM 

Relative 2 posit iolmenecic ltheplatiomn 
Input Z  positiongoteene 1th plattoun 


Lnpue 2 yOS 1c POmwOt mec nls 
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DATA SIMULATION PROGRAM DATA DECK SETUP 


ie 


Wrewcatamaceh=10, Use of the JSP 1s segregated into 


four basic parts. 


Card 


a. (Three cards) IFM Platform Information 
baa eeOne «care } Emitter Platfiorm Information 
Gee i Omea ic S ) Emitter Information 


(ome tive emitters) 


b. (One card) Emitter Piatctorm: Intormacien 


eu Clwo cards) Emitter Information 


(Up to ten platforms) 


ads (One card) End of data card with any negative 


number in columns 21-30 


The format for each required card is shown below. 


eet tiieelattorm Lntermat on 


Columns Format Variable 
1-10 110 NUNI 
11-20 IG NGAUS 

21-30 FLO Vie 
31-40 alors @ TO 
41-50 F10.0 TF 
51-60 F10.0 TINT 
61-70 F10.0 PHS 


oi 


Units 

Odd intcecer 

Reso Wit) Wetec 
Odd Integer 
Right Justified 
10° TDS/SEe 
Seeonads 

seconds 


Microsecconds 


dbm 





js Oue 9, 


Pang’. 0 


AL(0) 0 


EPOn0 


F10.0 
Jeidh(8) 18 
Faro" 0 
110 


XO 
YO 
ZO 
CSEO 


SPDO 
FMAX 
FMIN 
NBIN 


EM Cer. batt orm .htorlatron 


Columns 


41-50 


Sle sal 


Eimeeer intormacion 


Columns 


Formac 


Iie (8, 
Pale? 
Loz. 
|i 3 8, 


EeOr 0 
i 


Format 


lad). 10 
11118) 50, 
FEO. 0 


Ren a 


Vie ata ee 


XI 
YI 
tos 


oat 


SPpa 


NMTR(TI) 


Variable 


FREQ 
PRI 

PW 
PTRANS 
BEAM 


Nautical miles 
Nautical miles 
Fee 

Degrees clockwise 
fvom- North 

Knots 

MHZ 

MHZ 

ice cen 


Rigme CuStitied 


Units 

Nautical miles 
Nautical miles 
reece 

Degrees clockwise 
from North 

Knots 


Prvtte @ ean 


Unrate 

MHZ 
milliseconds 
microseconds 
Kilowatts 


Degrees 





Sl 00 GEOG, GAIN db 


61-70 Fia0:. 0 SCAN Seconds per revo- 
lution 
Z 1 jie IF(I) JLeejziers 
Z tel Penh} Integer 
3 Tl IPW(I) ineeces 
4 Ie ETE) Integer 
5 Ei Lr SCAniGn } Integer 


fmerna or data card 


Columns 21-30 Any negative number 


See hee program listed herein docs not include coding for 
bse Ot the parameter variations. This feature can easily 
be added by using utility subroutines GAUS and RANDU or 
any other random number generation scheme. 

eimesprorgram listed herein contains coding for only two 
types Of antenna scan--circular and steady. Coding for 
additional types was being prepared when it was determined 
what the run times for the ATAC Simulations were going to 
Be. 2s the simulations take very long periods to run on 
modest data sets, the addition of the additional features 
was delayed. Coding to add additional scan types can be 
mead iyeaccomplished. 

Peet neetinits of the input variables were selected so that 
TomeseGould be coded without the nmeéd for usang the 
eapemential format. It is felt that this feature will help 


meduee data deck format errors by allowing the input data 


0m 





to be placed anywhere within the 10 column field as long 


as a decimal place is used somewhere. 
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